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The Low Background infrared Calibra- 
tion Facility (LBTR) at the Narional In- 
stitute of Standards and Technology 
has been in operation for calibration 
measurements of ihc radiant power 
emitted from infrared radiation (IR) 
sources, such as cryugenic blackbodies, 
for more than 2 yeat^. The IR sources 
are sent to NIST by customers from in- 
dust^y'j govemmcni, and university labo- 
ratories. An absolute ciyngenic 
radiometer is used as the standard dc- 
lector to measure the total radtcint 
power at its aperture. The low back- 
ground is provided by a closed cyde 
helium refrigeration system that mam- 
tains the inner parts of ihc calibration 



chamber at 20 K. The radiance temper* 
ature of the blackbody is deduced from 
the measured power and compared 
■mXh the btackbody temperature sensor 
data. The calibration procedures and 
data analysis are illustrated using Ihc 
measurements of a typical blackbody. 
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1. Introduction 

The national need to establish a primary stan- 
dard reference for calibrating infrared sensors led 
to the establishment, in 1989, of the Low Back- 
ground Infrared Calibration Facility (LBIR) at the 
National Institute of Standards and Technology, 
This facility has been serving a variety of users from 
industry and government laboratories by providing 
calibrations of cryogenic blackbody sources. These 
sources are used by the customer for calibrations of 
their infrared sensor systems. An absolute cryo- 
genk radiometer (ACR) of Ihe electrical substitu* 
tion type has been developed as the standard 
reference detector for the LBIR calibrations. The 
ACR measures the total radiant power incident on 
its precision aperture. Its operating range for radi* 
ant power measuremenis is 20 nW to 100 ]xW with 



an expanded uncertainly of ^ 1% (95% confidence 
level). The physical description of the ACR and its 
characterization as an absolute detector for radiant 
power measurements have been described in Ref. 
[1]. The LBIR chamber is shown in Fig, L 

The cylindrical vacuum tank housing the ACR 
detector and the customer blackbody source is 
made of stainless steel and is 152 cm long and 60 cm 
in diameter. Two shields made of copper (shown in 
the cutaway portion of Fig* 1) separated by a 2.54 
cm gap are cooled by a helium refrigerator which 
circulates cooled helium gas (15 K) through copper 
lines vacuum brazed to these shields. The inner 
shield operates at 20 K and provides the low ther- 
mal radiation background inside the chamber. The 
inner surface of the inner shield has a highly 
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Fig. 1» LBIR chamber with partial cutaway showing the major features of the apparatus. 



absorbent black coating of IITRI Nffl-2200 paint' 
to reduce scattered radiation. The ACR can be 
placed in any of the three ports of the LBIR cham- 
ber, which are located at distances of approxJ- 
mately 30 cm, 65 cm, and 100 cm from the 
biackbody source under test. The placement of the 
ACR can be chosen so that the radiant power at its 
aperture is in the ACR measurement range. A 
complete description of the facility can be found in 
a previous publication [3]. 

The general requirements for biackbody calibra- 
tions are described in Sec. 2 of this article. The 
procedures for handling the cryogenics of the cali- 
bration setup and the procedures adopted for mea- 
surements and analysis are described in Sec. 3, A 
recent biackbody calibration illustrates the proce- 
dures. The biackbody in question has two sensors 
and is referred to as BB in this paper. 

2* Calibration Requirements 

Calibration requirements fall into two 
categories: 1) needs established by the customer 

and 2) needs of the NIST LBIR facility. A typical 



* Certain commencjal equipment, insfruments, or materials are 
KJentified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and 
Technobgy, nor ^oes it imply that the materials or equipment 
identified are ne<^ssari!y the best available for the purpose. J 



customer requirement for the BB is as follows: 1) 
measure the radiant power at the ACR aperture 
for the biackbody temperature range of 200 K to 
400 K in increments of 25 K for the source aper- 
ture diameter of 650 ^im and 2) calibrate the tem- 
perature sensors of the source by comparing with 
the deduced radiance temperature. NIST issues 
guidelines to the customer for preparing the black- 
body source to conform to the vacuum and allow- 
able contaminant level specifications of the LBIR 
chamber. The NIST requirement specifies that the 
partial pressure of hydrocarbons (beyond approxi- 
mately 45 u) should not exceed L33 x 10"^ Pa for a 
total base pressure of 133x10"^ Pa. Other re- 
quirements are described in detail in Ref, [3], 
NIST also supplies the customer with a mounting 
plate and guidelines for installing the plate on the 
biackbody source [3]. 

3. Calibration Procedures 

3*1 Cryogenic Apparatus and Preparation for 
Measurements 

The biackbody is integrated into the LBIR 
chamber after testing its vacuum integrity. Cryo- 
conditioning of the LBIR chamber follows with the 
inner shield of the chamber being cooled to ap- 
proximately 20 K and the outer shield being cooled 
to approximately 40 K- The biackbody mounting 
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shelf and the main chamber have separate cold he- 
lium feeds with flow control valves to allow diver- 
sion of additional coolant to the blackbody, if 
needed. The cryoshields have several silicon diode 
thermometers attached at various points. These 
thermometers are monitored during calibration to 
ensure stability of the IR environment. 

In order to provide a steady background flux for 
the ACR, a copper disk separates the blackbody 
chamber &om the ACR chamber. This disk^ called 
the isothermal plate, is moderately thermally insu- 
lated from the inner shield. The isothermal plate is 
fltted with an aperture that blocks the background 
radiation from the blackbody chamber without lim- 
iting the field of view of the ACR aperture. Com- 
plete isolation of the ACR chamber is achieved^ if 
necessary, by a shutter in front of the isothermal 
plate that can be operated by remote controL The 
shutter is thermally anchored to the isothermal 
plate with a flexible copper braid. The temperature 
of the isothermal plate and shutter are monitored 
by three precision silicon diode thermometers* The 
isothermal plate is controlled to 0*05 K by using an 
integral heating element and the temperature is 
held at 22 K. 

3.2 Measurements 

Figure 2 shows the calibration setup inside the 
LBIR chamber, Blackbody sensor data and radiant 
power data are collected at each blackbody tem- 
perature and aperture setting specified by the cus- 
tomer. Data are collected and recorded at 
approximately one second intervals for a period of 
three minutes at each setting. These data form the 
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Fig. 2. Blackbody calibration setup mstde the LBIR chamber. 



basis for estimating the repeatability of the instru- 
mentation. Cryogenic preparation of the ACR and 
the complete calibration sequence are repeated on 
three different days to test the reproducibility of 
the entire measurement system. 

3.2.1 Blackbody Temperature Sensors The 
blackbody temperature is generally measured by 
platinum resistance thermometers (PRTs) 
mounted on the blackbody core, Voitages across 
the PRTs measured by a multimeter, and currents 
supplied by a constant current source are con- 
verted to resistance values using Ohm's law. The 
resistances are converted to temperature values if 
the customer provides data on the calibration of 
the PRT sensors. Nominal settings on the external 
temperature control electronics of the BB and 
measured temperatures for the BB sensors are 
shown in coSumns 1, 2, and 4 of Table 1, Experi- 
mental standard deviations from the 1 s repetitions 
at each setting on each day are shown in columns 3 
and 5 of the same table, 

3.2.2 Radiant Power The radiant power at 
the ACR aperture is measured as follows; The ra- 
diation from the blackbody is blocked from the 
ACR receiver by using the shutoff position on the 
blackbody aperture wheel (if available). The elec- 
trical heater power to the ACR receiver is set to a 
value higher (preferably 20% to 30%) than the ex- 
pected radiant power from the blackbody so that 
the temperature controller servo (ac bridge) can 
maintain a constant temperature at the ACR cone 
by varying the amount of electrical power. Prelimi- 
nary estimates of radiant power at the ACR aper- 
ture are made using the Stefan-Boltzmann law and 
the approximate distance between the ACR aper- 
ture and the blackbody aperture. After the ac 
bridge is balanced, the ACR is allowed to receive 
the blackbody radiation by turning the aperture 
wheel to a chosen aperture position. If the cus- 
tomer blackbody is equipped with a fixed aperture, 
the 152 mm (6 in) shutter in front of the isothermal 
plate is used for blocking the blackbody radiation. 
The difference between the initial electrical power 
setting with the shutter closed and the final electri- 
cal power setting with the shutter open is the radi- 
ant power at the ACR aperture. Data is taken 
continuously at one second intervals for at least 3 
min with the shutter open, in order to determine 
the repeatability of the instrumentation, 

3*2*3 Geometry The measurement of the ra- 
dius of the blackbody aperture, n (see Fig. 2)^ cor- 
rected for thermal contraction at 20 K, is supplied 
by the customer* The radius of the ACR aperture, 
rj, has been measured by the Precision Engineering 
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Table 1. 


Experimental data 














Nomfnal 


Sensor 1 


Std, dev. 


Sensor 2 


Std. dev. 


Measured 


Average 


Std. dcv. of 


(K) 


(K) 


(K) 


(K) 


(K) 


power (nW) 


power (nW) 


average (nW) 


200 


199.S97 


0.002 


199.281 


0.002 


71.9 


71.5 


0,5 (0.6%) 




199.S74 


0.005 


199.332 


0.005 


72.0 








199.947 


0.002 


199J28 


0.003 


70.6 






225 


224,780 


0.003 


224.266 


0.004 


113.5 


1133 


0.1 (0.1%) 




224.728 


0.003 


224.276 


0.006 


113.4 








224,S23 


0.CK)4 


224.306 


0.004 


113.0 






250 


249.662 


0.004 


249.268 


0.004 


1713 


172.0 


0.4 (0.2%) 




249,639 


0.007 


249.288 


0.008 


172.6 








249.733 


0.005 


249.322 


0.004 


172.2 






275 


274.685 


0.003 


274.408 


0.002 


252.6 


252.5 


0.2(0.1%) 




274.644 


0,006 


274.394 


0.007 


252.2 








274.719 


0.006 


274.443 


0.004 


252.8 






300 


299 J31 


0.007 


299329 


0.007 


35S.4 


358.0 


0.4 (0.1%) 




299321 


0.003 


299.346 


0,004 


3573 








299388 


0.005 


299.3S9 


0.004 


358.4 






325 


324.412 


0,005 


324346 


Q.(m 


497.6 


494.6 


1.4 (0.3%) 




324.457 


0.005 


324.375 


0.005 


492.8 








324.446 


0.003 


324379 


0.004 


493.4 






350 


349.306 


0.003 


349.378 


0.003 


665.8 


664.8 


0.8(0.1%) 




349337 


0.003 


349.395 


0.011 


663.7 








349.445 


0.006 


349.479 


0.006 


664.8 






375 


374.047 


0.0O5 


374.284 


0.005 


882.2 


878.3 


2.0(0.2%) 




374,069 


0.003 


374.292 


0.003 


875.7 








374.171 


0.006 


374.377 


0.007 


877.0 






400 


399.026 


0.005 


399.362 


0,005 


1138.0 


1137.1 


13(0.1%) 




399.056 


0.004 


399,379 


0.004 


1134.6 








399.131 


0.004 


399.444 


0,003 


1138.8 







Division at NIST at room temperature in three dif- 
ferent orientations. The average radius is corrected 
for thermal contraction as explained in Ref. [1]. 

The distance^ S, between tiie blackbody aperture 
and the ACR limiting (precision) aperture is mea- 
sured in two parts. The user measures the distance 
between the blackbody aperture and the front sur- 
face of a reference tab on the blackbody mounting 
plate as required by NiST [3], and NIST personnel 
measure the distance from the front surface of the 
tab to the ACR aperture. The chamber measure- 
ment is assisted by a Kaman proximity sensor [4] 
located internally to determine the final location of 
the reference tab to the ACR aperture before the 
chamber is evacuated. The measurements are 
made at ambient temperature and corrected for 
thermal contraction due to cooling of the chamber 
and the ACR to 20 K and 2 K, respectively. The 
radius of the aperture in the isothermal plate, r^. 



shown in Fig, 2, is chosen to limit the background 
flux from the blackbody front surface. However, it 
will allow the full cone of light from the blackbody 
aperture to reach the ACR aperture. The geomet- 
rical measurements for the BB are given in Table 
2. 

3,2,4 Radiance Temperature The following 
equation deduced from the Stefan-Boltzmann law 
is used to convert the radiant power to radiance 
temperature, T: 



■=k 






1/4 



where 



A I crini 






(1) 
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Table 2, Geometric data for BB calibration 




Meas. 
value 






Std. dev. 
mean 


L Radius (r^) of the ACR aperture at 2.2 K 
Z Distance between the apertures (5) 
3. BB aperture radius (n): 


1.4971 cm 
30.77 cm 






5x10-^ cm 
0.042 cm 


Measured at 
Nominal ambient 
mm mm 


Dcdnccd 
20K 
mm 


Std. 


dev. mean 
mm 


0.325 0.3251 


0.3244 




0.0013 


t 





and^i = -7^^ is the area of the hiackbody aperture; 
P is the radiant power and om is the Stefan-Boltz- 
mann constant, and equal to 5.67051 x 10"* W m"^ 
K"*. The expression for the configuration factor, 
Fj, given above is taken from Ref. [5] and to a first 
order approximation reduces to tttiVS^. The quan- 
tities n, Kit and S have already been defined, and 
their measured values at room temperature are ad- 
justed to cryogenic temperatures using the stan- 
dard reference data for contraction of materials 
[6]. 

3-3 Analysis 

3.3,1 Diffraction Correction It has been es- 
tablished in the literature [7, 8] that the radiation 
from the blackbody incident at the ACR aperture 
will not be solely determined according to geomet- 
rical optics because of diffraction effects at both 
limiting and nonlimiting apertures in the beam 
path. Therefore, diffraction losses at each one of 
the apertures in the beam path are estimated by 
usmg the procedures published in Refs, [7] and [8], 
The measured radiant power values and the de- 
duced radiance temperatures are corrected for the 
total diffraction Joss which is a product of the frac- 
tional Josses at each aperture in the beam path. 
The calculated correction factors (AP/P) given as a 
percentage of the measured power values for the 
BB are listed in Table 3 for various temperatures. 
The corresponding correction factor (AT/T) for 
the radiance temperature is approximately 1/4 of 
the correction (AP/P) for the radiant powen 



Tabl« 3, Calculated diffraction correction {AP) given as a per- 
centage of tbe measured radiant power {F). Systematic uncer- 
tamty (Type B) in (be correction due to approximations in the 
calculations '\^ + 10% of AP/P. 



Nomina I 


Diff, corr. 


Uncer. (1 <r) 


blackbody 


AP/P 




temp. { K ) 






200 


13% 


±0.13% 


225 


1.2% 


±0.12% 


250 


1.1% 


±0.11% 


275 


1.0% 


±0.1% 


300 


0.9% 


±0.09% 


325 


0.8% 


±0.08% 


350 


0.8% 


±O.0S% 


375 


0.7% 


±0.07% 


^POO 


0.7% 


±0.07% 



33-2 Measurement Uncertainties The uncer- 
tainties in the measurements are analyzed accord- 
ing to guidelines under development by the 
International Organization for Standardization 
(ISO) [9]. Uncertainties are caused by either 
systematic or random effects. Uncertainties from 
systematic effects (here after called systematic un- 
certainties for simplicity) are of the following types 
for this experiment: 1) uncertainties in theoretical 
corrections applied to measured data such as the 
uncertainties associated with the diffraction correc- 
tions; 2) uncertainties in the measured values of 
parameters that remain constant during the 
calibration. Examples of such uncertainties are the 
parameters that characterize the ACR response [1] 
which do cot change during the ACR power 
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measurements. These uncertainties are Type B un- 
certainties according to Ref [9] Subclause 333. 
Uncertainties from random effects (hereafter 
called random uncertainties for simplicity) arise 
from the statistical variation in the measurements 
during the calibration experiment. This type of 
variation is inherent in the blackbody PRT sensor 
measurements of temperature and the ACR mea- 
surements of radiant power during the calibration 
experiment. These are Type A uncertainties [9]. 

Uncertainty in Radiant Power 

The uncertainty in the ACR measurement of the 
radiant power is analyzed as follows: For each 
nominal temperature, radiant power data (Pi) are 
collected every 1-1/2 s for approximately 3 min, 
leading to approximately 120 data points. An aver- 
age, P, and standard deviation, Srpt, are calculated 
from this data. The standard deviation, Jxpt, is a 
measure of the repeatability of the measurements 
in each set, and in general, the ACR data show vety 
good repeatability. The measurement process is re- 
peated three different times to estimate a standard 
deviation, ^,pr, for reproducibility. The average over 
the three runs is calculated as; 

J-l 

a standard deviation (s,) is calculated from 

3 



sf=2 {Pi-pfn. 

/-I 



(4) 



The standard deviation {st) is a measure of both re- 
peatability and reproducibility of measurements 
[10] according to 



J?=5^+5^^/120. 



(5) 



Although more repetitions of the entire measure- 
ment process would be desirable for better statis- 
tics, time constraints limit the number of 
repetitions. The radiant powers^Pj, and the average 
powers, P^ for the BB are given in columns 6 and 7^ 
respectively, of Table L The standard deviations of 
the means, sJ^/X are given in column 8- 

Two sources of systematic uncertainty are evalu- 
ated as Type B following Ref. [9]. The standard un- 
certainty (i.e., 1 o- level) in the diffraction correction 
to the measured radiant power is based on diffrac- 
tion theory used to calculate the correction. This 
standard uncertainty is estimated to be ± 10% [7, 8] 
of the calculated correction based on scientific 



judgment (Sec. 4.6 Ref. [II]), as shown in column 3 
of Table 3 for the BB. The standard uncertainty in 
the characterization of the radiometer as given in 
ReL [1] is 0,12% of the measured radiant power. 
TTie standard uncertainty due to systematic effects 
is then the square root of the sum of the squares of 
the two components. 

The random and systematic uncertainties are 
listed separately in the calibration report. The com- 
bined standard uncertainty (Uc) is reported as the 
square root of the sum of the squares of all of the 
components. However, it should be noted that 
other methods of combining components of uncer- 
tainty from systematic and random effects are in 
practice for calculating the total uncertainty [12]. 
Therefore, all components of uncertainty are listed 
to give the customer the choice of using an alterna- 
tive method. 

The radiant power data from Table 1 after cor- 
rection for diffraction loss are shown in column 6 of 
Table 4. The uncertainty components associated 
with the average power are shown in columns 7 and 
8f and the combined standard uncertainty, Ue, is 
given in column 9. 

Uncertainty in the Blackbody Sensor Data 

The Type A (random) uncertainty in the black- 
body PRT sensor data is analyzed as follows: The 
PRT sensor data are collected in conjunction with 
the ACR power measurements. In general, the 
standard deviations for the PRT sensor data are ho- 
mogeneous and negligibly small* Repeated sets of 
data are averaged to obtain the mean and the corre- 
sponding standard deviation analogously to the 
method discussed earlier for the ACR power mea- 
surements. The mean temperature values of the BB 
with associated standard deviations are shown in 
columns 2 through 5 of Table 4. 

Radiance Temperature and Calibration Uncertainty 

The final radiance temperatures for the calibra- 
tion report are obtained by the least-squares analy- 
sis of temperatures, T^; (i = 1, .. .,ni J -1^ • * •, 3), 
deduced from Eq. (1), as a function of blackbody 
sensor readings, Xi^ (r ~1,. . .,n; / = 1,. _ . ,3), for 
each aperture setting. 

In order to evaluate a confidence band for the 
variability of the calibration curve, the following 
statistical procedure is adopted. The calibration 
equation is assumed to follow the model 

Tij=aQl-aiXii+a2Xi} ^a:^X,^j -h .,.. + 

a,Xtj-^€,^, (6) 
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TaMe 4. 
B8 


Average bJadcbody sensor data 


and rBdJani 


power corrected for diffract km 


1 effects 


at the ACR 


aperture for the 


Nominal 


Sensor 1 
data 


Srd. dcv. 
(K) 


Sensor 2 
data 
(K) 


Std, dev, 

(K) 


Average 
power 
nW 


Sjsl. 

uncer. 

nW 


Random 

uncer 

nW 


Comh. 

unccr. 

«c(l^) 

nW 


200 


199,906 


0.003 


199 J13 


0.003 


72.4 


0.1 


03 


03 


(0.6%) 


225 


224,777 


0.003 


224.282 


0.004 


114.6 


0.2 


OJ 


0.2 (0.2%) 


250 


249,678 


0.005 


249.293 


0.005 


173.9 


03 


0.4 


oi 


(03%) 


275 


274.683 


0.005 


274.415 


0.005 


255.0 


0.4 


0.2 


03 


(0.2%) 


300 


299347 


0.005 


299 J55 


0.005 


361.2 


03 


0.4 


0.6 


(0.2%) 


325 


324.429 


0.005 


324367 


0.005 


49S,6 


0.7 


L4 


1.6 


(0.3%) 


350 


349362 


0.006 


349.429 


0.007 


669.8 


1.0 


0.8 


13 


(0.2%) 


375 


374,120 


0,005 


374331 


0.005 


8843 


1.2 


2.0 


23 


(03%) 


400 


399.071 


0-005 


399.395 


0.004 


1144.6 


1.6 


13 


2.1 


(0.2%) 



where ao,* - * ^ dk are to be estimated, and £,), the 
random errors associated wiih the measurements 
Tij, are assumed to be independent with heteroge- 
neous variances u^j. The a,j are estimated by Si/ 
from approximately 120 data points for each nin. 
Weighted least squares [13] acxounts for the het- 
erogeneity of variances where the weights, >Vi, are 
calculated from the empirical variances by 



Wij = l/Sx? . 



(?) 



The standard deviations associated with the Tij 
are an order of magnitude larger than the standard 
deviations associated with the Xi} data, and, there- 
fore» random uncertainties associated with the 
measurements Xij are assumed to be negligible for 
the purpose of the least squares analysis. The de- 
gree, k^ for the polynomial in Eq. (6) is determined 
by a goodness-of-fit test [14] which compares the 
agreement among groups of three runs with the 
overall fit to the data. The lowest degree polyno- 
mial which satisfies the goodness-of-fit criterion is 
taken as the calibration curve. 

Given a future blackbody sensor setting, Xh, its 
calibrated radiance temperature value is given by 

n=ah^aiXf.-^a'2Xi^a',X^'h .,.,^-a^XL (8) 

where ^o, * • » * dk are least-squares estimates from 
calibration data. The random component of uncer- 



tainty associated with the predicied value is com- 
puted as 






(9) 



The constant F{95;k-\-l, n -k~l)is the upper 95 
percentile of Snedecor's /"-distribution with k +1 
degrees of freedom in the numerator and n—k-\ 
degrees of freedom in the denominator, and s (7\) 
is the standard deviation of the predicted value^ Ts , 
This uncertainty, based on the Working-Hoetelling 
confidence bands for the calibration curve, is valid 
for all future applications of the calibration curve 
as long as the model holds [15], The details of the 
statistical analj'Sis can be found in Ref. [15, 16]. A 
software package, OMNITAB [17], is used for the 
statistical computations. 

Data for the least-squares analysis of the BB are 
shown in Table 5. Temperature as measured by 
PRT sensor 1 is used as the independent variabfc, 
X. The analysis confirms that a linear function 
{k - 1) is sufficient for describing the data, and the 
following equulion 



n^ai+a\{XH) 



(10) 



gives the predicted radiance temperature (7i) for a 
blackbody temperature {Xk) measured by PRT 
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Table 5. 


Data for the kast-squares analysis of the BB 










Radiance 


Standard 




Nominal 


Senior 1 


temperature 


deviation 


Weight 


(K) 


(K) 


(K) 


(K) 






X 


T 


s 


W 


200 


199.897 


20K427 


3J0O 


0.0798 




199.874' 


201.539* 


1.420 


0.4957 




199.947 


200.546 


2.130 


02201 


225 


224.780 


225.736 


2.990 


0.1116 




324.728 


225.666 


1.384 


0.5224 




224.823 


225.5(KS 


1.975 


0.2564 


250 


249.6(52 


250.135 


2.100 


0.2278 




249.639 


250,593 


1J66 


0.736 




249.733 


250,463 


1.865 


0,2875 


275 


274.6849 


275,575 


1.102 


0.8239 




274,644 


275,466 


lJi9 


0.5751 




274.719 


275.605 


1.101 


0.8244 


300 


299J309 


300.683 


0.715 


1.9545 




299J211 


300.462 


hl35 


0.7760 




299,588 


300.683 


0.887 


1.2702 


325 


324.412' 


326.342- 


0.830 


1.4529 




324.457 


325.600 


0.603 


2.7502 




324.446 


325.658 


h076 


0.8641 


350 


349.306 


350.931 


0.644 


2.4G84 




349337 


350,657 


0.997 


1.0056 




349.445 


350.807 


0.869 


1.3233 


375 


374.047- 


376.500' 


0.660 


2,2930 




374.069 


375.766 


0.816 


1.5003 




374.171 


375.915 


1.407 


0.4647 


400 


399.026 


401.171 


O.S72 


1.3164 




399.056 


40(1.866 


1.266 


0.6242 




399,131 


401.240 


1.217 


0.6748 



" Data with standardized residuals from the Jinear fit exceeding 2.5. These outliers are noi included in 
the final analysis. 



sensor 1. The estimated coefficients and associated 
standard deviations are: 

fl5=L0O7;5{flJ) = 0001 . 

The percentage uncertainty in the deduced radi- 
ance temperature^ ST/T, is given by Eq. (1) and the 
theory of uncertainty propagation [18]. An approxi- 
mation on a Taylor series expression gives the rela- 
tionship between the variables as 

hTIT - hFf4P + bnOri + Sr2/2^2 + 55/25 . (U ) 



The first term on the right hand side of Eq. (U) 
shows the relationship lo the contribution from the 
uncertainty in the radiant power which has both 
Type A and Type B components as discussed above. 
The last three terms in Eq. (U) show the relation- 
ship to the contributions from the uncertainties in 
the measurement of the geometrical quantities, n, 
r:, and 5. A summary of the percentage uncertain- 
ties in the geometrical measurements (Tabic 2) and 
their propagated uncertainties for radiance temper- 
ature measurements are given in Table 6. 

The predicted NIST radiance temperature from 
the calibration equation for each average value of 
temperature measured by PRT sensor 1 is given in 
Table 7 for the BB calibration. The uncertainly 
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TabEc 6, Uncertainties in ihc mciiSJirenients of geometrical quantities for the BB 



Quantity 



Uncertamty In 
measurement (1 cr) 



Propagait^d unc^r. 
for temperature 
measurement (! er) 



1) Distance between apertures (BSfS) 0,135% 

2) Radius of ACR aperture (Br2/r^) 0.003% 

3) Radius of BB apenurc (^ifri) 02% 



0.06a% 
0.002% 
0.1% 

Total* 0.12% 



■ Square root of sum of squares of 1, 2, und 3 given above. 



T«b*e 7. Predicted radiance temperatures and uncertainties for BB 



Nominal 


Sensor 1 


Predicted 


ACR 


Systematic uncertainties 


1 (T random 


E3tpandcd 




data 


radiance 
temp 


char. 


Geometiy 
meas. 


Diff. 
cal. 


Total 
b 


unccrtainlJes 


uncer. (U) 




(K) 


(K) 


CK) 


(K) 


(K) 


(K) 


(K) 


(K) 


200 


199.906 


200J 


0.06 


0.24 


0.07 


0.25 


0.10 


0.7 (0.4%) 


250 


224.777 


225.4 


0.07 


0.26 


0.07 


7§ 


o.os 


0.7 (0J%) 


250 


249-67S 


250,4 


0.08 


0.30 


O.OS 


0.32 


0.07 


0.7 (0.3%) 


275 


274.683 


275.6 


0.09 


0,33 


0.05 


0.35 


0,05 


0.8 (0J%) 


300 


299.547 


300.7 


0.09 


0.36 


0.06 


OJS 


0,04 


0,8 (0.3%) 


325 


324.429 


325.S 


0.10 


0.39 


0.07 


0.41 


0.04 


0.8 (0.3%) 


350 


349.362 


350.9 


0.11 


0.42 


0.07 


0.44 


0.04 


0.9 (0.3%) 


375 


374.120 


375.8 


0.11 


0.45 


0.08 


0.47 


0.05 


1.0 (0J%) 


400 


399.071 


400.9 


0.12 


0.45 


0.08 


0.5 


0.07 


LO (0J%) 



components due to systematic effects are listed 
separately based on Eq. (11): 1) the uncertainly in 
the power SP/4Py due to the characterization of 
the ACR as an absolute detector [1] is listed in 
column 4; 2) the uncertainty in defining the geome- 
try as given in Table 6 is shown in column 5; and 3) 
the uncertainty in diffraction calcuJations due to 
approximations is shown in column 6. The total sys- 
tematic component, b, which is the root sum of 
squares, is given in coJumn 7, The Type A compo- 
nent, s(Th), of the uncertainty is obtained from 
least-squares analysis, and the expanded uncer- 
tainty, t/, as shown in column 9, is obtained by ex- 
panding Eq. (9) as follo\\'s [19]: 

(/ = [(Jt +1) F(95; fc + 1, fi -A --1)]'^ 



The multiplying factor in Eq. (12) for/i =27, k-\ 
for the BB experimental data is 2*6, The value 
shown in parenthesis in column 9 is the expanded 
uncertainly given as a percentage of the measured 

temperature* 

The statistictil procedure that leads to Eq. (12) 
did not take into account possible correlations be- 
tween Type B uncertainties because of the relation- 
ship between variables shown in Eq, (11), However, 
calculations of expanded uncertainty, (/, using the 
complete covariance matrix and least square fitting 
using appropriate weights have been carried out 
and the results are found to be the same as given by 
Eq, (12). The calibration constants in Eq. (10) also 
are found to be essentially same as reported earlier. 

Figure 3 shows the temperature value from PRT 
sensor I plotted on the^-axis and the measured ra- 
diance temperatures (Table 5) and corresponding 
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predicted radiance temperature (Table 7) on the 
y-axis. The solid line in Fig* 3 connects the pre- 
dicted radiance temperature values and represents 
Eq, (10), Figure 4 shows the 95% confidence band 
for the difference between the calibrated tempera- 
ture and the PRT temperature. The expanded un- 
certainty given in Table 7 is within 0.5%, 
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Fig. 3* Measured temperatures and calibrated radiance tem- 
perature from polynomial fit plotted as a fuiK;tk>n of blackbody 
temperature setting using PRT sensor 1. 
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250 275 300 325 350 
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F%. A. Solid [me: DtSerence between the deduced radJance 
temperature and the temperature measured by the blackbody 
PRT sensor 1 plotted as a function of blackbody temperature 
setting. The points represent the same for measured radiance 
temperatures and are not all of equal weights. Dashed 
lines; Upper and lower bounds of 95% confidence bands. 

However, it is always an open question why for 
example, at the nominal blackbody temperature of 
the 400 K does the temperature sensor 1 measure 
399.071 K whereas the radiance temperature is 



400,9 K as shown in Table 7* The calibration curve 
does not address the effects of systematic errors 
such as unaccounted scattered light entering the 
ACR, aperture warming of the blackbody etc. How- 
ever, experimental designs to eliminate such sys- 
tematic effects are planned and implemented 
within the time allotted for customer calibrations. 
The experimental checks show these effects should 
be negligible. For future calibrations, spectral in- 
strumentation is also being added to address this 
question further by measuring the emittance from 
the blackbody as a function of wavelength between 
2 |xm to 30 jjjn and comparing it with the currently 
assumed value of unity* 



4< Summary 

The LBIR facility is in operation for customer 
blackbody calibrations using the ACR as the abso- 
lute detector. The experimental procedure and 
analysis of data that generate the calibration report 
for the customer are described in detail, and a re- 
cent blackbody calibration at the facility serves as 
an example of the procedure. Calibration uncer- 
tainties in radiance temperatures of less than 1% 
have been realized for typical flux levels of 25 nW 
and above at the ACR aperture. Experimental pro- 
cedures to reduce systematic uncertainties are con- 
tinuously investigated and implemented to improve 
the accuracy of measurements. Future improve- 
ments of noise control in the ACR operation could 
allow measurements down to 10 nW at the ACR 
aperture with 1% uncertainty. Capability for mea- 
suring spectral emittance of customer sources in the 
2 to 30 micrometer wavelength region is being de- 
veloped. 
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Disclaimer 

References made in this paper to particular 
brand names or specific suppliers of a service are 
made for ease of understanding by the reader and 
do not constitute an endorsement of products or 
senice by the National Institute of Standards and 
Technology over other competitive suppliers of 
similar products or service. 
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